Suppression of Magnetic State Decoherence Using Ultrafast Optical Pulses 
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It is shown that the magnetic state decoherence produced by collisions in a thermal vapor can 
be suppressed by the application of a train of ultrafast optical pulses. 



In a beautiful experiment, Itano et al. demonstrated 
the Quantum Zeno effect Q . A radio frequency pi pulse 
having a duration on the order of 250 ms was applied to 
a ground state hyperfme transition. At the same time, 
a series of radiation pulses was used to drive a strongly 
coupled ground to excited state uv transition. The r/and 
strong transitions shared the same ground state level. 
Itano et al. showed that excitation of the rf transition 
could be suppressed by the uv pulses. They interpreted 
the result in terms of collapse of the wave function - 
spontaneous emission from the excited state during the 
uv pulses is a signature that the uv pulse projected the 
atom into its ground state; the lack of such spontaneous 
emission implies projection into the final state of the rf 
transition. This paper triggered a great deal of discus- 
sion, especially with regards to the interpretation of the 
results ||. 

A necessary condition for a quantum Zeno effect is a 
perturbation of a state amplitude on a time scale that 
is short compared with the correlation time of the process 
inducing the transition. In the experiment of Itano et al, 
this time scale is simply the duration of the pi pulse, 256 
ms. On the other hand, if one wished to inhibit particle 
decay or spontaneous emission j3| , it would be necessary 
to apply perturbations on a time scale that is short com- 
pared with the correlation time of the vacuum, an all 
but impossible task. In this paper, we consider the inhi- 
bition of collisional, magnetic state decoherence, by the 
application of a train of ultrafast, optical pulses. This 
correlation time of the collisional perturbations resulting 
in magnetic state decoherence is of order of the dura- 
tion of a collision and is intermediate between that for 
the coherent pi pulse applied by Itano et al. and the 
almost instantaneous, quantum jump-like process pro- 
duced by the vacuum field. It should be noted that re- 
lated schemes have been proposed for inhibiting decoher- 
ence in systems involving quantum computation Q, but 
the spirit of these proposals differs markedly from the 
one presented herein. 

The rapid perturbations of the system are a necessary, 
but not sufficient, condition for a mechanism to qualify as 
a Quantum Zeno effect. The perturbations must involve 
some "measurement" on the system for the "Quantum 
Zeno" label to apply. The suppression of magnetic state 
coherence discussed in this paper is not a Quantum Zeno 
effect in this sense. We will return to this point below. 



We envision an experiment in which " active atoms" in 
a thermal vapor undergo collisions with a bath of foreign 
gas perturbers. A possible level scheme for the active 
atoms is depicted in Fig. [I]. At some initial time, an 
ultrashort pulse excites an atom from its ground state, 
having angular momentum J = 0, to the m = sub- 
level of an excited state having J = 1. The duration of 
the excitation pulse r p is much shorter than the dura- 
tion of a collision r c (r c is typically of order 1 ps). As 
a result of elastic collisions with the ground state per- 
turbers, population in the J = 1 sublevels equilibrate at 
a rate r co ; that is typically of order 10 7 — 10 8 s^ 1 per 
Torr of perturber pressure. The transfer to the m = 1 
substate is probed by a circularly polarized pulse act- 
ing on the J = l,m = 0— ► J = excited state transi- 
tion, applied at a time T~ o j following the initial excitation 
pulse. For the sake of dcfinitcness, we assume that the 
perturber pressure is such that equilibration occurs in a 
time r~; =0.1 — 1.0 ns. The question that we address 
in this paper is " How can one inhibit this magnetic state 
decoherence by subjecting the active atoms to additional 
external radiation fields?" 
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FIG. 1. Energy level diagram. The collisional interaction 
couples the magnetic sublevels in the J = 1 state. 
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As was mentioned above, the key to any Zeno-type ef- 
fect is to disrupt the coherent evolution a system from its 
initial to final state. In our case, the coherent evolution 
from the initial m = states to the final m = ±1 states 
is driven by the collisional interaction. Thus it is neces- 
sary to disturb the system on a time scale that is short 
compared with the collision duration r c . To do this, we 
apply a continuous train of ultrashort pulses that cou- 
ple the m = level to the excited state having J = 
shown in Fig. |l|. The pulses are assumed to have dura- 
tion t p <C t c and are assumed to be off-resonance; that is 
the atom-field detuning is large compared with t" 1 . As 
such, each pulse simply produces an ac Stark shift of the 
m = sublevel of the J = 1 state, resulting in a phase 
shift of this state amplitude. As a consequence, the exter- 
nal pulses break the collision-induced, coherent evolution 
of the atom from its initial m = state to the m = ±1 
states. If the pulse strengths are chosen such that the 
phase shift is a random number, modulo 2ir, and if many 
pulses occur during the collision duration t c , then the 
atom will be frozen in its initial state. It is interesting 
to note that collisions, which are normally viewed as a 
decohering process, must be viewed as a coherent driving 
mechanism on the time scales considered in this work. 

To obtain a qualitative understanding of this effect, it 
is sufficient to consider a model, two-level system, con- 
sisting of an initial state |0) (corresponding to the J = 1, 
m = state) and a final state |1) (corresponding to the 
,7=1, m = 1 state, for example). The Hamiltonian for 
this two-state system is taken as 

h = v c (t) (|o) (i| + |i) <o|) + nJ2 A.(tiM(* - *0 |o> (o| 

i 

(i) 

where V c (t) is a collisional perturbation that couples the 
two, degenerate states, and A s (ij) is the ac Stark shift 
of state |0) produced by the external pulse occurring at 
t = t L . For simplicity, we take V c (t) to be a square 
pulse, V c (t,b) = h(3(b), for 0< t < r c . The quantity b 
is the impact parameter of the collision. Without loss 
of generality, we can take the collision to start at t = 0. 
The collision duration t c can be written in terms of the 
impact parameter b characterizing the collision and the 
relative active atom-perturber speed u as T c (b) = b/u. 
Moreover, to simulate a van der Waals interaction, we 
set (3(b) — (C/6q) (°o/b) 6 , where C and b are constants 
chosen such that 2C/(6qu) = 1. The quantity bo is an 
effective Weisskopf radius for this problem. An average 
over b will be taken to calculate the transition rate. 

The external pulse train is modeled in two ways. In 
model A, the pulses occur at random times with some 
average separation T between the pulses. In model B, 
the pulses are evenly spaced with separation T. In both 
models, the pulse areas A s (^)t p are taken to be random 
numbers between and 2tt. A quantity of importance is 
the average number of pulses, no = r c (bo)/T = bo/(uT), 
for a collision having impact parameter bo- 



A. Randomly-spaced pulses 

The randomly spaced, radiative pulses act on this two- 
level system in a manner analogous to the way collisions 
modify atomic electronic-state coherence. In other word, 
the pulses do not affect the state populations, but do 
modify the coherence between the levels. The pulses can 
be treated in an impact approximation, such that dur- 
ing a collision, the time rate of change of density matrix 
elements resulting from the pulses is poo = pu = and 



Pio/pio = Poi/poi 



■r 1 



-iA s (U)T f 



(2) 



where T = T _1 is the average pulse rate and we have 
used the fact that the pulse area is a random num- 
ber between and 2tt. Taking into account the colli- 
sional coupling V c (t,b) between the levels, one obtains 
evolution equations for components of the Bloch vector 
w = pu- poo = 2/9n -l,v = i(pi - poi) as 

dw I dx = U (y)v; dv / dx = —U (y)w — n(y)v , (3) 

where x — t/r c (b) is a dimensionless time, y = b/bo 



is a relative impact parameter, and U(y) 



— „-5 



and 



n(y) — noy are dimensionless frequencies. These equa- 
tions are solved subject to the initial condition w(0) = 
— 1; v(0) = 0, to obtain the value pu(x = no) = 
[w(x — 1, y) + 1]/2. The relative transition rate S is given 
by 



S(n ) = 2nNubl 



ydypu(x = l,y,n a )/2, (4) 



where N is the perturber density. A coefficient, R(no), 
which measures the suppression of decoherence, can be 
defined as 



R(no) 



ydypn(x = l,y,n Q )/ / y dy pn{x = 1, y, 0) 



Solving Eqs. (0), one finds 



pu(x = l,y,n Q ) = 



1 - 



r 2 - n 



e 1 e 2 



r-2 



ri,a = {-noy ± VM 2 - ^y- 10 ) /2. 



(5) 

/2; (6a) 
(6b) 



It is now an easy matter to numerically integrate Eqs. (|^) 
to obtain R(no). Before presenting the numerical results, 
we can look at some limiting cases which provide insight 
into the physical origin of the suppression of decoherence. 

A plot of pu(x = 1, y, no) as a function of y = b/bo is 
shown in Fig. ||for several values of uq. With decreasing 
y, pu increases monotonically to some maximum value 
Pu(y m ) and then begins to oscillate about pn = 1/2 
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with increasing amplitude. One concludes from such 
plots that two effects contribute to the suppression of 
coherence. The first effect, important for large no, is a 
reduction in the value of y m . The second effect, impor- 
tant for no of order unity, is a decrease in the value of 
pu(y m )' Let us examine these two effects separately. 




1.0 
y=b/b 

FIG. 2. Graph of p u as a function of y — b/bg for several 
values of no. For values < y < 0.45 not shown on the graph, 
pu oscillates about an average value of 1/2. For no ^ 0, the 
oscillation amplitude increases with decreasing y. 

For very large no, n^ 66 3> 1, one can approximate pu 
over the range of y contributing significantly to the in- 
tegral (§ as p n (x = l,y,n ) ~ (l - e-*~"/«°) /2. By 

evaluating the integrals in (^|), one finds a suppression of 
decoherence ratio given by 



R(n ) = 0.95/ng 



2/11 



(7) 



The n 2 ^ n dependence is a general result for a colli- 
sional interaction that varies as the interatomic sepa- 
ration to the minus 6th power. It can be understood 
rather easily. The pulses break up the collision into 
noy segments, each having a (dimensionless) time du- 
ration Xb = X/(noy)- Each segment provides a perturba- 
tive contribution to pu of order y~ 10 (n y)~ 2 , provided 
y < y w , where y w is to be determined below. The to- 
tal population from the entire collision interval varies as 
Pu ~ y~ 10 (noyy 2 n y = y~ lx jn a . Of course, p n can- 
not exceed unity. One can define an effective relative 
Weisskopf radius, y w , as one for which pu = 1, namely 
Vw = bw/bo = n Q The total transition rate varies as 

Vw ~ n Q 2 ^ 11 j m agreement with (Q). As no ~ oo, the 
atom is frozen in its initial state. 

For values of no of order unity, the dominant cause of 
the suppression of decoherence is a decrease in the value 
of pii(y m ), rather than the relatively small decrease in 
y m from its value when no = 0. For values no < 3, 
approximately 45% of the contribution to the transition 
rate 5(no) originates from y > y m , and, for these val- 
ues of n , y m ~ it' 1 ' 5 and Pll (y m ) ~ (1 + e -"«/ 2 ^ /5 )/2. 



This allows us to estimate the suppression of decoherence 
ratio as i?(n ) = [0.55 + .45(1 + e -"°/ 27rl/5 )/2], such that 
R(l) = 0.93, R(2) = 0.88, R(3) = 0.84. These values 
are approximately 70% of the corresponding numerical 
results, indicating that the decrease in pn(y m ) accounts 
for approximately 70% of the suppression at low no, with 
the remaining 30% coming from a decrease in y m . The 
first few collisions are relatively efficient in suppressing 
decoherence. With increasing no, the suppression process 
slows, varying as n 2 / 11 , In Fig. ^, the suppression of de- 
coherence ratio R(hq), obtained by a numerical solution 
of Eq. (|5j), is plotted as a function of no- 
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FIG. 3. Graph of the suppression of decoherence ratio R as 
a function of no for randomly and uniformly spaced pulses. 



B. Uniformly Spaced Pulses 

We consider now the case of equally spaced pulses, hav- 
ing effective pulse areas that are randomly chosen, mod- 
ulo 27T. The time between pulses is T, and no = T c (bo)/T. 
For a relative impact parameter y = b/bo, with m < 
n(y) = noy < m + 1, where m is a positive integer or 
zero, exactly m or m + 1 pulses occur. The effect of 
the pulses is calculated easily using the Bloch vector. At 
x = 0, w = — 1 and v = 0. The Bloch vector then un- 
dergoes free evolution at frequency U(y) = y~ 5 up until 
the (dimensionless) time of the first pulse, x s = t s /T c (b). 
The pulse randomizes the phase of the Bloch vector, so 
that the average Bloch vector following the pulse is pro- 
jected onto the w axis. From x = x s to x s +T/r c (b) = x s 
+l/n(y), the Bloch vector again precesses freely and ac- 
quires a phase UT — y~ 5 /n(y) — y~ e /no, at which time 
the next pulse projects it back onto the w axis. Taking 
into account the periods of free precession and projec- 
tion, and averaging over the time x s at which the first 
pulse occurs, one finds 

w(y) = [1 - n(y)] cos[?T 5 ] 

+n(y) / dx s cos[y~ 5 a; s ] cos[y~ 5 (l - x a )]\ 
Jo 

0<y< 1/no, 
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w(y) = [m + 1 - n(y)][(m + l)/n(y) - 

x / efe s cos[y~ 5 a; s ] cos m_1 [y~ 6 /n ] 

«/ 1— rri/n(y) 

x cos[y~ 5 {l - .t s - (m - l)/n(j/)}] 
+ Ky) - m] [1 - m/n(j/)] _1 

~l—m/n{y) 

dx s cos[y~ 5 x s ] cos m [y~ 6 /"-o] 

x cos[y~ 5 {l - x s - m/n(y)}]; 
tn/riQ < y < (m + l)/no f° r to > 1. (8) 

In the limit that no 3> 1 , for all impact parameters that 
contribute significantly to the transition rate, approxi- 
mately n(y) pulses occur at relative impact parameter y, 
implying that w(y) ~ cos n ^ [y~ 5 /n(y)] and 



R(n ) 



(l-cos n °y[y- 6 /n }) 

(1 - cos[y- 5 ]) 
(1 - [1 - y- 12 /2nl] n °y) 



(9) 



(1 - cos[y 5 ]) 

1 _ e -y- 11 /2n X 



(1 - cos[y 5 ]} 



0.84 

2/11 ' 



which is the same functional dependence found for the 
randomly spaced pulses. Note that the form {1 — 
cos™^ [y~ 15 / n (y)]} is identical to that found in theories 
of the Zeno effect §. 

The suppression of decoherence ratio R(jiq), obtained 
from Eqs. (||) and (||) [using pn = (1 + w)/2], is plotted 
in Fig. 3. The fact that it lies below that for randomly 
spaced pulses is connected with the difference in the av- 
erage collisional phase shift acquired between radiation 
pulses for the two models. The oscillations in R(uq) ap- 
pear to be an artifact of our square pulse collision model. 
In the absence of the pulses, the first maximum in the 
transition cross section occurs for y max = (tt) -1 , corre- 
sponding to a 7r collision pulse. With increasing no, the 
pulses divide the collision duration into approximately 
n(y) equal intervals. If these pulse intervals are odd or 
even multiples of it, one can enhance or suppress the con- 
tribution to the transition rate at specific impact param- 
eters. Numerical calculations carried out for a smooth 
interatomic potential do not exhibit these oscillations. 



C. Discussion 

Although the collisional interaction has been modeled 
as a square pulse, the qualitative nature of the results 
is unchanged for a more realistic collisional interaction, 
including level shifts. In fact, for a smooth interatomic 
potential that allows for an increased number of radiation 
pulses over the duration of the collisional interaction, the 
suppression is slightly enhanced from the square pulse 
values. Although the pulses are assumed to drive only 



the J = 1, to = — ► J = 0, excited state transition, it is 
necessary only that the incident pulses produce different 
phase shifts on the J — 1, m — and J = 1, m = 1 state 
amplitudes. 

To observe the suppression of decoherence, one could 
use Yb as the active atom and Xe perturbers. The Weis- 
skopf radius for magnetic decoherence is about 1.0 nm 
H, yielding a decoherence rate of ~ 10 10 s _1 at 500 Torr 
of Xe pressure at 300°C, and a collision duration r c (6o) 
~ 2.5 ps. Thus, by choosing a pulse train having pulses 
of duration r p —100 is, separated by 0.5 ps, it is possible 
to have 5 pulses per collision. If an experiment is carried 
out with an overall time of 100 ps (time from initial exci- 
tation to probing of the final state), one needs a train of 
about 200 pulses. To achieve a phase shift A s r p of order 
2tt and maintain adiabaticity, one can take the detuning 
5 = 3xl0 13 s _1 and the Rabi frequency Q ~ lxlO 14 s _1 
on the J = l,m = — »J = 0, excited state transition [pj. 
The corresponding, power density is ~ 1.5 x 10 11 W/cirr , 
and the power per pulse is ~ 150 uJ (assuming a 1 mm 2 
focal spot size). This is a rather modest power require- 
ment. With 5 pulses/collision duration, one can expect 
a relative suppression of magnetic state decoherence of 
order 40%. 

Finally, we should like to comment on whether or not 
the effect described in this work constitutes a Quantum 
Zeno effect. Normally, the Quantum Zeno effect is pre- 
sented as a projection of a quantum system onto a given 
state as a result of a measurement on the system. In 
the experiment of Itano et al., this "measurement" is re- 
flected by the presence or absence of spontaneously emit- 
ted radiation during each uv "measurement" pulse. The 
measurement pulse must be sufficiently long to produce 
a high likelihood of spontaneous emission whenever the 
atom is "projected" into the initial state by the pulse. 
Following each measurement pulse, the off-diagonal den- 
sity matrix element for the two states of the rf transition 
goes to zero. In our experiment involving off-resonant 
pulses, the number of Rayleigh photons scattered from 
the J = level during each applied pulse is much less 
than unity. As such, there is no Quantum Zeno effect, 
even if suppression of magnetic state decoherence occurs. 
On average, each pulse having random area destroys the 
coherence between the J = 1, m = and J — 1, to = ±1 
state amplitudes, but does not kill this coherence for a 
single atom. With an increasing number of radiation 
pulses, no, however, both the average value and the stan- 
dard deviation of the transition probability tends to zero 
as Uq 1 for each atom in the ensemble. 

The experiment of Itano et al. could be modified to 
allow for a comparison with the theory presented herein, 
and to observe the transition into the Quantum Zeno 
regime. If the pulses that drive the strong transition are 
replaced by a sequence of off-resonant pulses, each pulse 
having a duration t p much less than the time, T w , re- 
quired for the pi pulse to drive the weak transition, and 
each pulse having an effective area, A s t p = (ri 2 /4(5)r p , 
that is random in the domain [0,27r], then the pulses will 
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suppress the excitation of the weak transition (it is as- 
sumed that Cl/S <C 1). If the upper state decay rate is 7, 
then the average number of Ravleigh photons scattered 
during each pulse is n = (fi/45) -fr p . For n < 1, there is 
suppression of the transition rate as in our case, while, for 
n > 1, there is suppression and a Quantum Zeno effect. 
There is no average over impact parameter, since exactly 
[Tn/T] or ([T ff /T] + 1) pulses in each interval between the 
pulses, where [x] indicates the integer part of x. 
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